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I. INTRODUCTION 



Heat conduction in one-dimensional systems has evoked a considerable amount of interest 
in recent years. In particular, the validity or otherwise of the Fourier law of heat conduction 
in such systems is a non-trivial question, and hence one that has been the subject of lively 
discussion A comprehensive account of the current status of the problem is provided in 
Ref. 0. 

Most of the results known in this regard are based on numerical studies of lattice models. 
An alternative and useful line of development is the analytical study of a model system, 
albeit a simplified one, which we call the Jepsen gas: a system of identical classical 
point-particles of mass m moving on a line and undergoing perfectly elastic collisions when 
neighboring particles meet^. This system is a special case of the more general model 
investigated in the context of the so-called adiabatic piston problem^], in which a central 
heavy particle of mass M is in a gas of particles of mass m on its left and right. The 
problem then is to analyze the dynamics of the central particle in the thermodynamic limit 
in which the one-dimensional gases to its left and right are in thermal equilibrium at specified 
temperatures and densities. The Jepsen gas is a singular limiting case of the unequal mass 
situation. The latter is not integrable, in marked contrast to the case M = m, which is 
integrable. Notwithstanding this simplification, non-trivial irreversible behavior is observed 
in the Jepsen gas in the limit N oo when appropriate averages over initial conditions are 
performed. In the initial studies of this systemj^, 0], several quantities of interest such as 
the statistics of the displacement and velocity of one of the particles, which we shall refer 
to as the central or tagged particle (or "piston", as it is the counterpart of the adiabatic 
piston in the model at hand), have been calculated exactly by performing an average over 
equilibrium initial conditions for the rest of the particles. In particular, it can be shown 
that the motion of this tagged particle becomes diffusive asymptotically, i.e., converges to 
Brownian motion. 

Recently, the model has been revisited and studied in greater detaii0,Q], one of the mo- 
tivating factors being its relationship to the adiabatic piston problem. The non-equilibrium 
situation implied by different velocity distributions for the gases to the left and right of the 
tagged particle has been analyzed. In general, this particle acquires, in the thermodynamic 
limit, a systematic drift velocity over and above its diffusive motion. A notable feature is 
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that this drift velocity is exclusively fluctuation-induced (6|. In the special case when the 
velocity distributions of the gases on the two sides are Maxwellian, so that the pressure of 
each gas can be identifled with ks times the product of its (linear) number density and 
temperature, an interesting feature emerges: even when the pressures of the two gases are 
equal, the drift velocity of the tagged particle does not vanish. Rather, it is directed from 
the lower temperature (higher density) side to the higher temperature (lower density) side. 
The drift velocity vanishes when the product of the number density and the square root of 
the temperature is the same on the two sides of the tagged particle, a condition already 
recognized in Ref. j^. 

In the context of heat conduction in one dimension, the question that arises naturally 
is whether the heat flux can be calculated for the Jepsen gas. We show in this paper that 
this problem can be solved analytically: it is possible to calculate exactly the amount of 
energy that is transfered through the tagged particle or piston. At flrst, it would appear that 
the issue is a trivial one in the same way as it is in a linear chain of harmonic oscillators. 
In the latter system, the energy just travels ballistically, being carried by phonon modes 
which do not interact with each other. In a similar fashion, the kinetic energy carried by 
any particle in the Jepsen gas is just transfered upon collision to the next particle, and 
hence moves ballistically along the line. As a result, a Fourier law, which predicts the 
diffusive spreading of thermal energy, does not appear to be valid. We will show by an 
exact and explicit calculation that this "hand- waving" argument, while formally correct, 
is nevertheless misleading in the sense that a non-trivial energy flux obtains in the model. 
A closed expression can be derived for this quantity, comprising two components. One of 
these involves the asymptotic drift velocity of the tagged particle, while the other is present 
even when the drift velocity vanishes. This permits the identiflcation, and hence a natural 
separation, of the "mechanical work" and "heat" contributions, respectively, to the energy 
transfer. In particular, for initial conditions corresponding to thermal equilibrium of the 
gases to the left and right of the tagged particle (with densities n~ and and temperatures 
T~ and T+, respectively) such that n~\fT~ = n~^VT+, so that the drift velocity is zero, the 
asymptotic steady-state energy flux through the tagged particle is shown to converge to a 
constant value. For small values of (T+ — T^) , this quantity is linearly proportional to the 
temperature difference itself, which can be interpreted as a manifestation of the Fourier law. 
Moreover, the "mechanical work" contribution itself is shown to have a part that is linear 
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in the drift velocity, as might be expected, as well as a nonlinear part that starts (for small 
drift velocities) with the fourth power of the drift velocity. 

We reiterate the following point. The Jepsen gas is admittedly a simplified special case 
of the more "realistic" models of one-dimensional transport that have been the subject of 
much attention. Nevertheless, it is the fact that analytical (and hence unambiguous) results 
can be obtained for this model that makes it worth studying, because these results help shed 
light on several of the essential issues involved in energy transport in one dimension. 

The plan of the rest of this paper is as follows. In the next section, we introduce the 
notation and summarize the salient results pertaining to the Jepsen gas for the purpose at 
hand. In Sec. IIIII fand the Appendix), an exact formula is derived for the energy transfer 
across the piston at any instant of time. This leads to a formula for the total energy flux in 
the stationary state, which is obtained and analyzed in Sec. IIVI In Sec. El the contribution 
to the energy flux coming from the stationary heat flux is identified, and the nature of the 
Fourier law for the model is clarified. The rate of mechanical work is also deduced, and its 
nonlinear dependence on the stationary drift velocity elucidated. Section ED contains a few 
concluding remarks. 



II. NOTATION AND RECAPITULATION 



It is helpful to recapitulate in brief the relevant features of the model, in the notation used 
in earlier work ^ Q]- The tagged particle, located at X = at t = with an initial velocity 
Vq, separates a gas of particles in the interval [—L, 0) on its left from a gas of A^"*" particles 
to its right, in (0, L]. Their initial positions Xj (where — < J < — 1 for the gas on the left, 
and 1 < j < A^^ for the gas on the right) are independently and uniformly distributed in the 
corresponding intervals. Their initial velocities Vj are drawn from normalized distributions 
4>~{V) and (f)^{V), respectively. To avoid unnecessary complications, we shall assume that 
these are symmetric distributions, i.e., 4>^{V) = (p'^l—V). It can be shown that the system 
has a thermodynamic limit in which A^+ ^ oo, L ^ oo with finite densities lim N^/L = n^, 
provided only that the mean speeds = dVV (p'^iV) are finite. Note that 0^ need 

not be Maxwellian distributions, although this is the case of direct interest in the present 
context of heat conduction. As all the particles (including the tagged particle) have equal 
masses, they merely exchange their identities on their original linear trajectories Xj{t) = 
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Xj + Vjt in each collision. This is what enables one to derive exact results for the system. 
Among other quantities, the exact one-particle distribution function of the tagged particle, 
P{X,V,t\0,Vo), can be found. The complicated stochastic motion of the tagged particle 
(henceforth termed the "piston") can be interpreted [3] as being driven by two independent 
Poisson processes (corresponding to collisions from its left and right, respectively) with state- 
and time- dependent intensities n~a{X/t) and (3{X/t), where 

POO pW 

a{W)= dV{V -W)(I)-{V), f3{W)= dV {W - V) (j)^ (V) . (1) 



JW J-oo 

The system does not equilibrate in the conventional sense of the term: the initial set of 
trajectories persists for all t. The piston acquires an asymptotic mean drift velocity W 
given by the unique solution of the implicit equation 1^ 

n-a{W) = n+(3{W) . (2) 

The asymptotic velocity distribution of the piston is a superposition of (p'^iV) and (p^{V), 
and is given by 

Pst{V) = [n-(P~{V)e{V -W)+ n+(P+{V)e{W - V)] I S(W) , (3) 

where 

;>oo pW 

S(W) = n- |a'(W)| +n+/?'(W) = n" /_ c/V0-(V) +n+ / c/V0+(\/) (4) 

JW J~oo 

is the normalization factor. Only in the special case of the homogeneous system, defined by 
n~ = = n and ip'^iV) = 4>~{V) = (piV), can one have stationarity and an approach to 
equilibrium, in the sense that W = and Pst{V) = 4'{V)- Note that W may vanish even 
in the inhomogeneous system, if the condition dV V (p^iV) = n+ dV V (p~^{V) (or 

'^"(l^l)^ — '^^(l^l)^) happens to be satisfied: that is, the mean rates at which the piston 
suffers collisions from its left and right, respectively, are equal. 

We point out at this juncture that the stationary distribution Pst{V) in Eq. (jS)) may in fact 
be anticipated on physical grounds, as follows. In the equal-mass case under consideration, 
we may regard the particles as simply "passing through" each other in each collision, because 
they are identical particles. In the stationary state, after the memory of the initial velocity 
of the piston has been lost, it will only encounter such particles via collisions from its left, 
as have velocities greater than its asymptotic mean drift velocity W*, whatever that be 



(trajectories corresponding to < W* from the gas on the left will never again intersect 
the piston's trajectory). Similarly, the piston will encounter only those particles from the 
gas on its right as have velocities less than W*. Moreover, after each collision, the piston 
simply acquires the velocity of the particle it collided with. Hence the stationary velocity 
distribution must necessarily be a superposition of the form 

Pst{V) oc n-(t)-{V) e{V - W*) + n+0+(V) e{W* - V) (5) 

apart from a normalization factor. The latter is easily seen to be just 1/2(1^*) where the 
function S is as given by Eq. Taking the first moment of this equation, we find that 

consistency demands that W* satisfy precisely the same implicit equation (Eq. as that 
written down above for the drift velocity W. In other words, not only is the form in Eq. 
(jni) for Pst{V) deducible on direct physical grounds, but also the existence of an asymptotic 
drift velocity W and the implicit equation satisfied by it. 

III. EXACT EXPRESSION FOR ENERGY TRANSFER 

To derive of an exact formula for the rate of transport of energy across the piston, we 
first calculate the difference AE{t) between the mean energy E{t) of the gas on the left of 
the piston at time t and its initial value E{0). 

We have 

^(0) = J:-Ln~ {V') = \^ T.^-N- {V' d{-X,)) , (6) 

while 

m = \m Ea E,^a {v,'0{Xa{t) - X,{t)) 5^^[iV-, Y.,^J{X,{t) - X.m) ■ (7) 

Here each sum runs over the entire set of particle labels {— A^~, ...,—1,0,1,..., A^"*"}; the 
subscript a is used to distinguish the instantaneous variables of the piston from those of the 
other particles; 5^''[n, m] is the Kronecker delta 5nm and (. . .) denotes an average over the 
already-specified distributions of the initial conditions comprising the set {Xj , Vj}. Using 
the crucial fact that 

Y,d{Xa{t)-Xe{t))=N- (8) 
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for all t > 0, the gain in energy of the gas on the left at time t can be written as 

= 5^5^(1^^ -X,(t)) 6''^[N', Y.^{X^{t)-X,{t))]). (9) 

The representation bnm = {2Tri)^^ dz enables us write this as 



AE{t) = 




where the contour encloses the origin. We break up AE{t) as 

AE{t) = AE\t) + AE-{t) + AE+{t) , (11) 

corresponding to the three distinct contributions to AE{t) coming, respectively, from the 
cases a = 0, — A^^ < ct < — 1 and 1 < a < A^^. The calculation of these quantities is broadly 
similar to, but a little more intricate than, that involved P, Q| in evaluating quantities like 
the one-particle distribution function P{X, V, t\XQ, Vq) of the piston. An outline of the main 
steps leading to the expression for AE{t) is given in the Appendix. The exact result for 
AE{t) is given by Eqs. (rO|l - (rO0|l . 



IV. FORMULA FOR THE STATIONARY ENERGY FLUX 

To find the asymptotic or steady-state rate of transport of energy across the piston, we 
pass to the long-time limit of the expression for AE{t). Using the asymptotic behavior 
Ir{z) ~ /{27rzY^'^ as \z\ — 00, it is evident from Eq. ()A.8|) that A£'°(t) decays exponen- 
tially to zero as t ^ 00, unless it so happens that n~ a{0) = n'^ (3(0), i.e., the asymptotic 
drift velocity of the piston is zero. In that case Ai?°(t) ~ t^/^ at long times; but this again 
implies a rate of transport that vanishes (like t~^/^) as t ^ 00. Therefore AE^(t) may be 
dropped from further consideration. 

Turning to Eqs. (|X9|l and (lA.lOj) for AE-{t) and AE+{t), we observe that the 
asymptotic behavior of the modified Bessel functions leads to the occurrence of a factor 
exp < —t (?2~a(w)) — (^n'^P{w)) > in the integration over w (which runs from —00 
to -|-oo). Since the exponent has a unique zero at precisely the asymptotic drift velocity 
W as defined in Eq. (0), the long-time behavior of AE^{t) can be deduced by a standard 
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Gaussian approximation about w = W. The leading asymptotic behavior of each of these 
two contributions is then seen to be ~ t, implying the existence of a finite, non-vanishing 
stationary rate of energy transfer lim^^oo dAE{t)/dt, that we denote by AEgt- Carrying out 
the algebra required, our final result for this quantity is remarkably simple (for reasons to 
be explained shortly). We find 

= im|n+ dU U\W - U) - dU U\U - W) 0-(f/)}. (12) 

For ready reference, we recall that the asymptotic drift velocity W of the piston is given by 
the implicit equation (j2|), which may be written in the alternative form 

_ ^ n- dU U 0- jU) + n+ JZ dU U 0+(^) ^^^^ 
n~J^ dU<p-{U)+n+jZ^ dU<PHU) 

We remark that the result in Eq. ()12|) is valid for arbitrary initial velocity distributions 
0^(1^), and not just for Maxwellian 0^ (see below). It must also be noted that the existence 
of the third moment of the speed, is necessary in order that AEst be finite. The 

exact result in Eq. (|T^ has several other noteworthy features, that we take up in turn. 

(i) Physical interpretation: The formula in Eq. ()12j) can be given a direct physical inter- 
pretation, extending the heuristic argument described in Sec. Ullfor PstiV) and W: once 
again, the fact that the particles merely exchange velocities upon colliding implies that the 
mean stationary rate of energy transfer across the piston can be constructed by an energy 
balance argument. If the piston has a stationary drift velocity W, the trajectories belonging 
originally to the gas on its right that collide with it are those corresponding to velocities 
from —oo to W. The energy carried by (a particle on) each such trajectory is ^mU^. The 
rate of collisions of such trajectories with that of the piston is {W — U), the second factor 
being the relative velocity between the two. Therefore 

+ rW 

^ / dUU^W -U)(j)+{U) (14) 

is the mean rate at which energy is transferred to the piston by the gas on its right, in the 
stationary state. Exactly the same argument shows that 
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dUU"^ {U -W)(t)-{U) (15) 
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is the mean rate at which the gas on the left transfers energy to the piston. The difference 
between the two is precisely the formula of Eq. (fT^ for the mean stationary rate of transfer 
of energy from the right to the left across the piston. Our rigorous derivation serves to 
corroborate this physical argument, in addition to providing an exact result for the time- 
dependent transients as well. 

(ii) Dichotomous and three-valued velocity distributions: The dichotomous velocity distribu- 
tion 

0+(V^) = <p-{V) = <PiV) = i [5{V + c) + 5{V - c)] (16) 

yields, as always, a simple and tractable special case that serves as a useful check on the 
calculations. We have, in this case, W = c{n~ — n~^)/{n^ + n^). Evaluating the various 
quantities appearing in Eq. (fT^ . we find that AEst vanishes identically for the dichotomous 
distribution above. However, valuable insight into the structure of the result for AE'st 
is provided by a slight generalization of the dichotomous distribution of Eq. |TB|) to the 
distribution 

= <p-{V) = <P{V) = ix5{V) + i(l - /x) [5{V + c) + 5{V - c)] , (17) 

where < fx < 1. It has been shown in Ref. 0] that the introduction of a non- vanishing 
probability mass fiatV = significantly alters the long-time properties of the homogeneous 
system. For instance, the velocity autocorrelation function acquires a t^^/^ tail, in contrast 
to its exponential decay in the case of a dichotomous (f){V). In the present context, too, 
a non- vanishing value of fi leads to a strikingly different result for AE^t- We have in this 
instance 

W = (18) 

(1 - fi)n> + (1 + /i)n< ' 

where n> = max (n^, n+) and n< = min (n~, ?2+). Computing the various quantities occur- 
ring in Eq. (fT^ . we arrive at the result 

2 (l-/i)n> + (l+^)n<' ^ ^ 

Thus AEst vanishes identically if n+ = = n (which, together with 0+ = 0~ as already 
imposed by Eq. p7|) . implies a homogeneous system), as it must in the homogeneous system. 
Similarly, AEst vanishes when /i = (the case of a dichotomous (p'^CV)), or when fi = 1 
(the trivial case of no motion at all). Comparing the expressions in Eqs. ()18|1 and |T9|l . we 
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observe that AEgt can be written in this case in the reveahng form 

AE,t = lmc'^fin-{-W). (20) 

(Recall that we have defined AE as the rate of transfer of energy to the gas on the left of 
the piston, and that W is negative if > n~ , i.e., if the piston drifts to the left.) 
{Hi) The case of Maxwellian distributions: The situation that is of direct physical interest 
is of course that of Maxwellian velocity distributions (p^iV) characterized by temperatures 
T^. For arbitrary densities and the asymptotic drift velocity W is still given by the 
solution of a transcendental equation, which reads in this case 

— n-(2fcBT-/m7r)V2e— ^'/2fes^- -n+(2fcBT+/m7r)i/2e-'"^'/2'=sr+ , , 

W = ^ ^ =2 • (21) 

n- [l - erf {mW /2kBT-y/^] + n+ [l + erf {mW /2fcBT+)V2] 
In terms of n"^, and W as given above, we get 

A^t = (27rm)-i/2|r2+(A;BT+)3/2e— -r2-(A;BT-)3/2e-'"^'/2'=^^"| 
+\W |n+(mW^ + kBT+) [l + erf {mW^ /2kBT+y/^] 

+ n-{mW^ + keT-) [l - erf {mW'^ /2kBT-f/^] | . (22) 

The nonlinear dependence of Ai^st on W is explicit in this formula. We shall comment 
further on this subsequently. It is also clear that AE^^ is not proportional to the temperature 
difference [T^ — T^), though one might perhaps naively expect such a proportionality based 
on an incorrect identification of AEst with the heat fiux, together with an application of the 
Fourier law to the system under discussion. 

(iv) Mechanical work and heat contributions: This brings us, finally, to a very important 
point. The general expression obtained in Eq. (fT^ for AEst incorporates the effects (on the 
energy transfer rate) of both the drift of the piston and its diffusive motion (or fluctuations 
about its mean position, namely, about Wt) . In broad terms, one might regard the respective 
contributions as the rate of mechanical work done upon the gas on the left of the piston [W 
being the rate of compression or expansion, depending on whether W < or W > 0), and 
the rate at which its entropy changes. However, these contributions are intricately mixed 
up in the formula for AEgf Moreover, the dependence of AEgt upon W is in general highly 
nonlinear. Disentangling these contributions will enable us, in principle, to isolate the parts 
of AEst that may be identifled with the "heat flux" AQst and the "rate of mechanical work" , 
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respectively. In some cases, the former may vanish altogether — as in the example of the 
distribution considered in Eq. ()17p . for which Eq. (j^UI) shows clearly that Aii^st arises entirely 
from the drift of the piston. In general, however, such a clear separation does not occur in 
the system under study. It is also important to bear in mind the fact that a non-vanishing 
drift velocity W is itself a consequence of the statistics of collisions in the system under 
consideration jjof. 



V. STATIONARY HEAT FLUX 



A direct way to isolate and examine the "heat flux" is to impose the condition of zero 
drift {W = 0) by adjusting, for instance, the value of the ratio n~ /n~^ of the densities of the 
two gases. As we shall now see, this leads to considerable simplification in the formula for 
AEst, which we are now justified in re-labeling as AQst- 

As emphasized more than once, the asymptotic drift velocity W of the piston vanishes if 
the densities and initial velocity distributions (p^ are related by the condition n~a{0) = 
n"'"/?(0), or n^{\V\)^ = n'^{\V\)~^, i.e., the mean rates at which the piston suffers collisions 
from the gases on either side of it are equal. The motion of the piston is then purely diffusive 
in the long-time limit, with a variance (A^(t)) that tends asymptotically to 2Dt, where the 
diffusion coefficient is given by 0, Q 

D = 2n-{\V\)-/{n- +n+)\ (23) 

We note that the homogeneous system (defined by n~ = n~'',0" = 0^) automatically has 
W = 0. The converse is not necessarily true, of course: W = does not necessarily imply 
a homogeneous system. Recall also that we have already considered situations in which 
(f)~^{V) = (f)~{V) = (piV), but n'^ n~; we then have an inhomogeneous system in which, 
moreover, W ^ 0. 

Setting W = in Eq. (fT2|l . we find that the stationary energy fiux in the absence of drift 
is simply 

A4t = lm[n-^{\V\Y-n-{\V\')-] , (24) 
in terms of the third moments 

POO 

{\V\Y = 2 / dUU^(l)^{U) (25) 
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of the particle speed in the two gases. For the homogeneous system, of course, AQst vanishes 
identically, as it must. 

Turning again to the case of Maxwellian distributions (p'^iV) at temperatures T^, the 
drift velocity W vanishes, as is well known, if n^\/T^ = n^VT^. Equation ()22j) then 
reduces to 

A4t = n+(^j A;b(T+-T-). (26) 

We see that the heat flux is indeed given by the Fourier law in this case, with a "coefficient 
of thermal conductivity" that is proportional to the square root of the temperature. 

It is revealing and instructive to pause at this stage to compare the result in Eq. ()2(j|l with 
that for the heat flux between two reservoirs held at different temperatures and densities, 
coupled by virtue of their sharing a common piston j^. When the mass of the latter is equal 
to that of the gas particles, the Boltzmann equation can be solved exactly for the stationary 
distribution; concomitantly, the heat flux conveyed via the shared piston from one reservoir 
to the other can also be calculated. While the exact expression for this quantity (see j^) 
is slightly more complicated than that of Eq. above, it reduces to the latter result to 
leading order in the temperature difference (T"*" — T~), apart from an extra numerical factor 
tt/^/2. This lends additional support to our identification of the portion of the energy flux 
that corresponds to the heat flux. The fact that the heat flux is about twice as large in 
the case of the shared piston is readily understood by recalling that the piston now suffers 
collisions with the gas particles belonging to both reservoirs, so that the effective rate of 
collisions is roughly twice as large. For a detailed discussion of the relevance of the problem 
of the shared piston to the general questions addressed here, we refer to P|. 

Finally, let us return to the full expression for the stationary energy flux AE^t in Eq. ()12p. 
and analyze its dependence on the drift velocity W. For sufficiently small W, we may expand 
AEst in powers of W, taking care to incorporate the fact that W = imposes the condition 
'^"(l^l)^ — ""-^(l^l)^ '^'^ the parameters occurring in the coefficients of the expansion. (This 
has been done in the Maxwellian case, Eq. (j2ni) above, to reduce n+(T+)3/2 - n-{T-f/^ to 
n"'"(T"'")^/^(T^ — T^).) For distributions (p^iV) that have derivatives of all orders at ^ = 0, 
the formal expansion in powers of W is given by 
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where 

^{V) = n+(f)+{V)~n~(j)-{V), (28) 

^C^) denotes its k^^ derivative, and AQst is given by Eq. (|24p . This representation isolates 
the contribution to AEst owing to the systematic drift of the piston from that arising from 
the fluctuations about its mean position. We observe that the part that is nonlinear in W 
is 0{W^). Considering the Maxwellian case once again, we find 

1/2 / r,_m+m- \ 1/2 



^st 



(^) n^ks{T^-T-)+(^-^f^) ks{n-Vf^-n^VT^) (29) 
\ ZTcm J \ Sirm j ^ ' 



correct to first order in the difference {rr\jT^ — W^^/T^^, the next term being of fourth 
order in this quantity. In contrast to the relatively simple expression to which AEst reduces 
when rT \JT^ = n^^/T^ (Eq. (j^Bj) ). no significant simplification of the general expression 
in Eq. (j^^ occurs when n~T~ = n'^T^ (equal pressures on either side of the piston). 
Some simplification does occur, however, when the temperatures on the two sides are equal, 
= = T (but n'^ ^ n"). We find that the drift velocity is now given by the implicit 
equation 

1/2 



W=( !}—^] ( ^] e-^w'/2ksT ^ {mW'/2kBTy/' 
\n~ + n+ J \ m-K J 



(30) 



The corresponding stationary energy flux is found to be 



\ Svr / \ m J 

Written in terms of the moments of the (Maxwellian) velocity distribution, this is just 

AEst = i m (TF - {V^)) in- - n+) {\V\) e--^^ I-^^^t _ (33) 

Finally, if the two densities are equal {n^ = = n, but 7^ T^), we find that W becomes 
independent of n, and AE^t is directly proportional to n. 



VI. CONCLUDING REMARKS 



We have already commented at the appropriate junctures on the special and interesting 
features of the structure of our analytical result for the stationary rate of energy transfer, as 
given by Eqs. (fT^ and ^T7\. We conclude with a comment on the question of the finiteness 
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or otherwise of the heat conductivity of the system under study [l|, It may be argued 
that the coefficient of heat conductivity of our system is essentially infinite: The contention 
is that in the expression for the heat fiux, the conductivity is the coefficient of the gradient 
of the temperature, and the latter is oc (T+ — T^)/L. This gradient tends to zero in the 
thermodynamic limit, although a finite energy fiux persists in the long time limit. Hence 
the coefficient multiplying the gradient has to diverge. 

While we agree with this argument, we point out that the system can be viewed in a 
different way: the gases on the left and right of the piston are heat reservoirs which have 
a single microscopic degree of contact, namely, the tagged particle we have termed the 
piston. There is no relevant length scale in the problem, so the "heat flux" is expected to 
appear solely as a result of this thermal contact. The fact that the Boltzmann calculation 
gives a comparable value for the conductivity supports the meaningfulness of our result. 
Furthermore, one can construct a physical system that will display exactly the behavior 
predicted by our model calculation. Consider a two- or three-dimensional cylinder of axial 
length 2L with a central piston that can move without friction along the axis of the cylinder. 
The compartments to its left and right are fllled with gases in equilibrium, at respective 
densities and temperatures n^,T^. In the limit of low densities and L oo, the dynamics 
of the piston will be exactly as described by the one-dimensional model: in the physical 
set-up, the particles only interact with the piston, while in our one- dimensional model, they 
exchange velocities upon collision, which is effectively tantamount to their passing through 
each other without interaction. 
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APPENDIX: CALCULATION OF AE{t) 



As stated in the main text, we write Eq. (fTUj) for AE{t) as the sum Ai?'^(t) + AE~{t) + 
AE^{t) of contributions coming, respectively, from the possibilities that at time t, the piston 
is (i) on its original trajectory, so that a = 0, or (ii) on a trajectory belonging to the gas on 
its left, so that —N~ < a < —1, or (iii) on a trajectory belonging to the gas on its right, so 
that 1 < a < A^+. Consider AE^ first. We have in this case Xq = 0, Vq = 0. We find 

AE\t) =ml^Y.{^' [^(-^^- - ^) - ^(-^^■)] > 

x\{[z + {l- z){e{-X, - V,t))] J] [1 + (z-i - l){e{-Xe - Vet))] (A.l) 

I- l'->r 

where the symbols I— and ^'+ are used to denote the fact that I and £' run over — < ^ < 
— 1 and 1 < < A^"^, respectively, in the products concerned. We shall use the convenient 
notation (• • • )^ for the corresponding averages. The evaluation of the contributions AE^{t) 
is somewhat more involved. We get 

^AT- (^|(Ar- - l){Vf [e{Xa + Vat-X,- V, t) - 9{~x^)])- 

+ ( [e{Xa + Vat-x^- V, t) - e{-x^)] )+} 

x[z + {l- z){9{X, + Vat-Xe- V,t));f'-' [z + {1 - z) 9{Xa + VJ)] 
x[l + {z-' - 1){9{X, + Vat-Xe- Vet)%f^y , (A.2) 

/ a 

while 

/I 
(^\^N-{V^^[9{Xa + Vat-X,- V,t) - 9{-X,)]Y^ 

+ (A^+ - l){Vf [9{X, + Vat-X^- V, t) - ^(-X,)])^} 

X [z + (1 - z){9{Xa + Vat-Xe- Vet));f' [1 + {z'' - 1) 9{Xa + K t)] 

X [1 + {z-' - l){9{Xa + Vat-Xe- Vet))^,f^-'y . (A.3) 

Evaluating the averages required, we get, for instance, 

{v^m-x, - V, t) - 9{-X,)]): = - /,7, dU r{u) - (t/L) dU r{u) , 

{V^ [9{-X, - V, t) - 9{-X,)]y = j:^/' dU <p^U) - {t/L) dU <p^U) , 
{9{-Xe - Vet)); = j'/'^dUriU) - {t/L) jy'dUUr{U) , 

{9{-X, - \/,,t)>; = j:^'dU<P^{U) - {t/L) jy^^dUU<j>^U) . (A.4) 
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The other averages are similarly calculated. Inserting all these results in Eqs. ()A.1|) - ()A.3|1 . 
we pass to the thermodynamic limit — > cxd, L ^ oo, such that lim N^/L = . Using 
the well-known expression for the generating function of the modified Bessel function Ir{z), 
the final expressions are as follows: 

Recall the definitions of the rates a{w) and P{w) in Eq. ((T]), and let 

\~{w ; n~) = rT a{w) , \^{w ; n^) = n'^ f^iw). (A. 5) 

Define the effective rates 

X{w ; n-,n+) = {X-X+y/^ , A{w ; n-,n+) = A" + A+ . (A.6) 

Further, let 

poo PW 

F{w]n-,n^)=n- I dU {w ~ U) (j)' (U) + n+ dU {w - U) (U) . (A.7) 

J w J — oo 

Then, suppressing the dependence in A and A for notational simplicity, we find 

AE\t) = imte-^(o)*/o(2A(0)t)F(0; (A.8) 

AE-{t) = ^mn^t^ f^^dwe-^'-'"^^\a'{w)\F{w; n-,n+) 

X {9{w) Io{2Xiw)t) + ei-w) [X+{w)/X-{w)Y/^ Ii{2X{w)t)} , (A.9) 

AE+{t) = ImnH'^ J^^dwe-^^""^' f3'{w)F{w] n-,n+) 

X {e{-w) Io{2X{w)t) +e{w) [X~{w)/X+{w)Y/^Ii{2X{w)t)} . (A.IO) 

AE{t) is the sum of the right-hand sides of Eqs. ()A.8jl - ()A.injl . 
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